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Scalability and coherence are two essential requirements for the experimental implementation of 
quantum information and quantum computing. Here, we report a breakthrough toward scalability: 
the simultaneous generation of a record 15 quadripartite entangled cluster states over 60 consecutive 
cavity modes (Qmodes), in the optical frequency comb of a single optical parametric oscillator. The 
amount of observed entanglement was constant over the 60 Qmodes, thereby proving the intrnisic 
scalability of this system. The number of observable Qmodes was restricted by technical limitations, 
and we conservatively estimate the actual number of similar clusters to be at least three times larger. 
This result paves the way to the realization of large entangled states for scalable quantum information 
and quantum computing. 
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Introduction The experimental implementation of 
quantum computing, driven by the promise of exponen- 
tial speedup for tasks such as the simulation of quan- 
tum physics [1] and integer factoring [ ] is a daunting 
challenge that requires exquisite levels of control over 
the quantum mechanical properties of numerous indi- 
vidual physical systems (quantum bits or, in this paper, 
quantum modes or Qmodes). The response to this chal- 
lenge spawned a wealth of experimental research efforts 
in widely different fields [ ] , striving to enable and main- 
tain quantum-coherent temporal evolution of quantum 
bits, while at the same time scaling up their number. 
Here, we demonstrate a breakthrough toward scalabil- 
ity: the novel, ultracompact implementation of quantum 
registers in the optical frequency comb (OFC) formed 
by the spectrum of a single optical parametric oscillator 
(OPO), thereby utilizing a capability for quantum infor- 
mation storage analogous to that exploited classically in 
FM radio or wavelength multiplexing. The classical OFC 
generated by ultrastable pulsed lasers has found ground- 
breaking uses in ultimate precision frequency measure- 
ments [4, 5]. In the case of the quantum OFC, each 
(Q)mode is well approximated by a quantum harmonic 
oscillator whose continuous-variable Hilbert space is de- 
fined by its amplitude- or phase-quadrature field observ- 
able (analogues of the position and momentum observ- 
ables). There is no known fundamental impossibility to 
the implementation of quantum computing with Qmodes 
[6-8], even though the implementation of quantum er- 
ror correction appears likely to require Hilbert-space dis- 
cretization [9, 10]. A method to create a frequency- 
degenerate A^-Qmode register was proposed, by use of 
N OPOs and a 2A'-port interferometer [ ], and demon- 
strated [12, 13] for 3 and 4 Qmodes. However, it was 
also shown that a square-grid continuous-variable clus- 
ter state of arbitrary size, suitable for universal one-way 



quantum computing [8, 14], can be generated in the OFC 
of a single OPO [15, 16]. In this work, we achieved the 
first step toward this goal: the parallel generation of 15 
quantum computing registers, each comprising 4 Qmodes 
in a quadripartite cluster state, in the quantum OFC of 
a single OPO. Requirements for the generation of larger 
entangled states include the experimental progress made 
in this work, along with a richer pump spectrum and a 
more tailored nonlinear interaction [ , ]. 

Experimental method The quantum OFC was gen- 
erated by a bowtie ring OPO containing two x-cut 
KTiOP04 (KTP) nonlinear crystals, of 10 mm length, 
and rotated by 90° from each other about the x axis. 
This ensured the perfect overlap of the respective OFCs 
of orthogonal linear polarizations y and z. One crystal 
was not phasematched. The other was periodically poled 
with two distinct periods: 9 jiia over a 3 mm length, and 
458 /im over 7 mm. The former quasiphasematched the 
zzz parametric downconversion, where the first letter de- 
notes the polarization of the pump field at frequency 2uJo 
and the other letters denote the polarization of the n*^ 
signal field pair at co'±n=^oi(^ — 1/2) A, with A=945. 66 
MHz the free spectral (FSR) range of the OPO cavity. 
The latter period quasiphasematched the yzy and yyz 
interactions simultaneously (dispersion was negligible for 
our values of n). The pump polarization was carefully 
adjusted in the {yz) plane, using OPO characterization 
by resonant second harmonic generation [^V], to yield the 
Hamiltonian [18] 

H = ihKY^ (^^-n,^ 4,^ + ^-n^y ^n^z + «-n,^ 4,?;) ^^.C. 
n 

(1) 

where ^ is the creation operator of the /c-polarized 
Qmode of frequency ujj . This Hamiltonian entangles the 
OFC as depicted in Fig.l and proven by the solutions of 
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FIG. 1. Principle of the experiment. The OFC of a single 
OPO is made polarization-degenerate by using a cavity with 
2 identical crystals rotated 90 degrees from each other in the 
polarization plane. One crystal simultaneously phasematches 
the zzz, yzy, and yyz nonlinear interactions, bottom. This 
creates square weighted cluster states, top, in blue [ ]. 

the Heisenberg equations for the n*^ Qmode quartet: 

Q+ = {[Q-n,y - Qn,y] + ^ [Q-n,z - Qn,.]} e""^ (2) 



P+ = {[P-n,y + Pn,y] + ^ [P-n,z + Pn,z]} - 



(3) 



Q- = {^[Q-n,y^Qn,y]-[Q-n,z^Qn,z]}e * (4) 
P_ = [P^n,y - Pn,y] - [P-n,z - Pn,z]} C ^ (5) 

where Q = a -\- and P = i{a^ — a) are the ampli- 
tude and phase quadratures, r is the squeezing parame- 
ter, and <l> = (>/5 + l)/2, which is the golden ratio. These 
four squeezed (quantum-noise reduced) field quadratures 
coincide, to local quadrature phase shifts left, with the 
nullifiers (entanglement witnesses) of a weighted square 
cluster state [18] (Fig.l) in the (unphysical) limit of infi- 
nite squeezing, where the cluster state is a zero-eigenvalue 
eigenstate of the nullifiers. The exponentiated nullifiers 
are thus the stabilizers of the entangled state [ ] (in 
contrast to the Qbit case, weighted Qmode cluster states 
are stabilizer states [ ]). For a pure state, observing the 
squeezing of a nullifier suffices to prove that the state 
has been prepared into a stabilizer state. For a sta- 
tistical mixture, the situation is more complicated but 
one can still use the van Loock-Furusawa criteria [ ] to 
prove quadripartite nonseparability. We experimentally 
demonstrated both. 

The setup is described in Fig. 2. The OPO, pumped at 
532 nm by a frequency- doubled, ultrastable continuous- 
wave NdiYAG laser (Innolight Diabolo), consisted in 
a cavity with low- loss mirrors and a 5% output cou- 
pler. The quantum OFC was separated into its orthog- 
onal polarizations, and quadrature combinations, e.g. 
Q-n,y^Qn,y Eqs. (2-5), wcrc measured by two-tone 
balanced homodyne detection with local oscillator (LO) 
fields at uj±n' The LO originated from another NdiYAG 
laser (Lightwave Electronics) which was phaselocked to 
the 1064 nm pump laser before it was frequency doubled 
[17]. The LO laser frequency could therefore coincide 
with ujo (Fig. 3), or differ from it for experimental veri- 
fications (Fig. 4). The uj±n frequencies were then gener- 
ated by phase electro-optic modulation (EOM) and sub- 
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FIG. 2. Experimental setup with HWP, half waveplate; 
QWP, quarter waveplate; FI, Faraday isolator; PZT, piezo- 
electric transducer; PLL, phase-lock loop; AOM, acousto- 
optic modulator; EOM, waveguide electro-optic modulator; 
KTP, KTiOP04; RTA, RbTiAs04 EOM; PBS, polarizing 
beam splitter; ATT, RF attenuator; SA, spectrum analyzer. 

sequent bandpass filtering by an optical cavity of same 
FSR as the OPO, in order to remove the carrier and sec- 
ond harmonics. The homodyne visibilities were 97% for 
the y polarization and 96% for z. Finally, the homo- 
dyne photocurrents from 95% efficient InGaAs photodi- 
odes (JDSU ETX500T) were preamplified and combined 
by RF splitters and attenuators in order to yield the vari- 
ance of nullifiers (2,5), observed synchronously on two 
spectrum analyzers while the LO optical path Olo was 
scanned. The measured observables can be expressed in 

■a^e'^^ as 



terms of generalized quadratures A{0) = ae 

A±{e)={[A.n,y{0)TAn,y{-e)] 

± [^_„,,(0)TAn,.(-^)]}e-'-*' 



(6) 



where phase values ^=0,7r/2 yield amplitude and phase 
quadratures, respectively. Note that the squeezing 
is independent of 0. Because we use two-tone ho- 
modyne detection, 6> is a function of 6lo and of 
the EOM phase Oo (Fig. 2), and the different nulli- 
fiers Eqs. (2-5) are obtained for the respective val- 
ues ((9LO,6>o) = (0,0);(0,7r/2);(7r/2,7r/2);(7r/2,0), modulo tt. 
Additional checks were made using LO polarization [17]. 

Results Fifteen sets of 4 Qmodes were measured for 
n=l to 15. The measurement results are displayed in 
Fig. 3. As can clearly be seen, the level of squeezing is 
constant over the whole set of 15 observed clusters, which 
establishes scalability in the OFC. Moreover, the maxi- 
mum value of n = 15 was not fixed by the quantum state 
preparation process but by measurement limitations: the 
14 GHz bandwidth limit of the EOM. The state prepa- 
ration bandwidth is given by the phasematching band- 
widths of the nonlinear interactions. We calculated these 
for the zzz and yzy processes, respectively, to be 616 and 



3 



Mill IIIIIIIII II JIIIIIIIII HIM 




llll l ll l ll l| lll | lll |l ll l lllllll 




M il l lll j ll ! lll|llllll | l l l m il 




0.0 2.5 



IIII I IIIIII M IIIIIIIIIII I IIIIII 





0.0 2.5 5.0 

Local Oscillator Phase (rad) 




0.0 2.5 5.0 



iii n i ii ii | iiii | ii M| i ii ii i iiii 




ll ll lll | III M II|l l lllll | l l l l lll 




0.0 2.5 



llll j l l llllllll|lllllll l ll | llll 




0.0 2.5 



I I IIIIIIIIIIIIIIIIIII I IIII I II J I 




0.0 2.5 5.0 

Local Oscillator Phase (rad) 




MM 11 Mill 



III I III M|I IIII|III II| I I IIIIIII 




I M III M I I I II II|IIIIIII IMI IIII 




-2.5 0.0 2.5 



lll | lll l lllllll|lll ll l ll lll | lll 





0.0 2.5 5.0 

Local Oscillator Phase (rad) 



FIG. 3. Scaling quadripartite entanglement in the op- 
tical frequency comb of a single OPO: nuUifier variance 
/\A±[e(eLo)Y, relative to the vacuum noise, versus the LO 
phase Olo- Note that the squeezing depends on Olo but not 
on ^ [ ]. Single-sweep measurements were taken at 1.25 MHz 
frequency; RBW: 30 kHz; VBW: 30 Hz. The Qmodes (black 
lines) measured are marked by the LO sidebands (purple ar- 
rows). The green line references half the pump frequency. 

47 GHz at 99% of the maximum. This indicates that a 
1% squeezing decrease occurs for the yzy interaction at 
47 GHz (a 10% decrease at 153 GHz). This is too weak 
an effect to be observed at our current squeezing level and 
cluster states n = 16, ...,47 should therefore be identical 
to the ones measured on Fig. 3. We therefore expect that 
three to ten times as many cluster states were generated 
than the 15 that were accessible with our setup. We can 
also flat-top shape the phasematching curve [ ] in order 
to optimize scalability. 

Phase-locking our two 1064 nm lasers to each other al- 
lowed us to make the crucial checks necessary to test the 
validity of our experimental results, in particular of our 
two-tone homodyne detection. These checks consisted of: 
using a single LO sideband, placing the LO sidebands 
at uncorrelated frequencies, and detuning the pump fre- 
quency from our quantum OFC. Figure 4 shows typical 
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FIG. 4. Detection and entanglement checks. As in Fig. 3, 
the plots display nullifier variance measurement relative to 
the vacuum noise versus the LO phase Olo- Top row, left 
and center, single-sideband detection displays no single-mode 
squeezing in the OPO comb, in any quadrature. Top row, left, 
and bottom row, left and center, LO sidebands coincide with 
uncorrelated comb lines, which yields no multimode squeezing 
whatsoever, again no matter the LO phase used. Bottom 
row, right, OPO pump detuned from the comb (2uOo 7^ 00 -n + 
ujn) which makes the nonlinear interaction singly (instead of 
doubly) resonant and yields negligible squeezing. Notations 
and legends as in Fig. 3. The dashed purple line references the 
phase-locked laser's frequency. 

results, which all agree with theoretical predictions [17] 
and clearly show no quantum correlations whatsoever, 
in stark contrast with the nullifier squeezing signals of 
Fig. 3. An essential point here is that all these verifica- 
tion results were insensitive to the LO phase, unlike the 
nullifier measurements. 

We finally address pure state preparation. The fact 
that the antisqueezing magnitude is larger than the 
squeezing one points to the existence of losses and ad- 
ditional classical noise (from the pump laser), and there- 
fore to the creation of a statistical mixture rather than a 
pure state. This can be alleviated by filtering the pump 
field with a "mode-cleaner" cavity, which we didn't do so 
as to maximize the pump power and hence the amount 
of squeezing. Nullifier squeezing is enough to claim en- 
tanglement if the state is pure. In order to ascertain this, 
we measured the squeezing spectra of A± at the optimum 
phases. Fig. 5. As can be seen, the state is pure for mea- 
surement frequencies above 5 MHz [ ], which validates 
our cluster-state preparation claim. 

In the case of a mixture, as in the case of our 1.25 
MHz measurement frequency (which yields more squeez- 
ing), one can use the van Loock-Furusawa (vLF) sep- 
arability criterion [6] in order to show that no Qmode 
can be placed in a factorized density operator of its own. 
A detailed analysis [ ] led to five vLF inequalities that 
must all be experimentally violated. While some of these 
inequalities have bounds at the vacuum level and are 
trivially violated by mere nonzero squeezing, others have 
bounds below the vacuum noise and therefore present 
higher violation thresholds. Figure 6 displays experimen- 
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FIG. 5. Squeezing spectra for left, and right. The 
squeezed trace on the left was recorded simultaneously with 
the antisqueezed trace on the right, for Olo = vr/2, and vice 
versa, for Olo =0, as with measurements in Figs. 3 & 4. 



tal results for the two most difficult such cases, which 
were clearly violated, thereby proving quadripartite en- 
tanglement even in a mixed state. 
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FIG. 6. Variances (AA+)^, blue, and (AA_)^, red, with mea- 
surement gains set to 0.3 in lieu of 1/$ = 0.618. (The vLF 
criterion allows us to deviate from nullifier measurements to 
find the optimum gain values for maximum violation.) The 
dashed lines indicate the sub-vacuum-noise violation levels 
required to prove inseparability. 

Conclusion We demonstrated that the optical fre- 
quency comb of a single optical parametric oscillator 
lives up to its promise as an extremely scalable sys- 
tem for quantum information. We simultaneously gen- 
erated a record number of quadripartite cluster states, 
in a record number of Qmodes, all equally entangled. 
The quantum comb was read by two-tone homodyne de- 
tection. Even though the size of the entangled states 
themselves is not a record, compared to the 14-ion GHZ 
state [23], we demonstrated stringent state preparation 
requirements for cluster states, a universal quantum com- 
puting resource. A practical quantum computer will re- 
quire an increase in both the number of entangled modes 
and amount of squeezing. However, the projective mea- 
surements required for one-way quantum computing can 
already be performed on the clusters that we gener- 
ated [24]. Variants of our setup will allow the genera- 
tion of multiple cube graphs [ ], and a scalable quan- 
tum wire and square-grid lattice [lo, i.]. We thank 
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EXPERIMENTAL GENERATION OF SQUARE CLUSTER STATES IN THE OPTICAL FREQUENCY 

COMB 

Hamiltonian 

Our optical parametric oscillator (OPO) contained two nonlinear crystals, placed in a two-waist ring cavity. A 10 
mm long periodically poled KTi0P04 (PPKTP) crystal was placed in the smallest waist. Three mm of the crystal 
were poled for quasiphasematching the zzz ("type-0") parametric downconversion, while the remaining 7 mm were 
poled for quasiphasematching the yyz/yzy (type-II) parametric downconversion, where the first letter denotes the 
polarization of the annihilated pump photon and the last two letters denote the polarization of the emitted signal 
(entangled) fields. The poling periods of 9 jam for zzz and 458 jam for yzy were chosen to ensure that both interactions 
displayed maximum efficiency at exactly the same temperature. These poling periods were calculated by using the 
Sellmeier equations given in Ref. [ ], with the temperature dependences from Ref. [ ]. We operated at 32.34°C, which 
was chosen by observing the temperature yielding the largest difference frequency interaction for both yzy and zzz. 
An unpoled x-cut KTP crystal, rotated by 90°, was placed in the other OPO waist. This nonphasematched crystal 
compensated for the birefringence of the PPKTP crystal, thereby making the free spectral range (FSR) of both y- 
and z-polarized light an identical A = 945.66 MHz and allowing for the spectral overlap of these polarization modes. 
Modes of both polarizations for each frequency pair were then linked together by parametric downconversion to form 
square cluster states. 

In the interaction picture, we can write the Hamiltonian of the triply concurrent process taking place in the OPO 

H = ifiKjyzy ^ttL^^^a^^^^ + o}j^ji^ya}_^^^ + ^ o}_^^^a}j^^^^ + i^.c, (7) 



where 



and 



^yzy 



2ttuj 

XyzyPy \P) 



^ ^zzz XzzzPz'^y dyzyLyzy^z'^y ^g'^ 

^yzy Xyzy (^y'^z d^zz^zzz (^y'^z 

assuming confocally focused Gaussian beams in the nonlinear crystal [ ]. The optical frequency modes {a-ri,yi ci-n,z) 
and {a-^n,y^CL-\-n,z) are positioned on opposite sides of the comb, with frequencies U-n = Uo — {n — 1/2)A , Un = 
co'o + (n — 1/2) A , respectively, where A = 945.66 MHz is the free spectral range of the OPO. 
We can write the Hamiltonian (7) in the general form 

H = ihhZyzy ^ Gij a\a}- + H.c, (10) 

where G is of Hankel form and called the adjacency matrix of an 1-L (Hamiltonian)- graph 

/O 1 \ 

10 

G- 1 e • ^^^^ 

Vi ^ 0/ 
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Solving the Heisenberg equations for Qmode evolution under we obtain the following squeezed joint operators, or 
nullifiers: 

{t)± 

where r = nyzyt and 

= (13) 

When ^ = 1, = is the golden ratio) and the real parts of these operators, phaseshifted by ^, coincide with 

the nullifiers A±{0) of the main text. 

The experimental setup is displayed in Fig. 2 of the main text. The pump polarization was set to an angle of 19.6° 
with respect to the y axis of the poled crystal, which meant that 89% of the light was y polarized and 11% was z 
polarized. This pump polarization orientation was used to counteract the difference in efficiency between the two 
processes, thereby setting ^ = 1, and 8± = It was determined that the SHG efficiency of the zzz interaction in 

our OPO was 4 times stronger than that of yzy by modeling a scan of resonant second harmonic generation (SHG) 
output vs infrared (IR) input polarization. Therefore, in order to match the interaction strengths of the zzz, yzy, 
and yyz interactions, it was necessary to pump the cavity with an 8:1 ratio of y light to z light [^]. 



Phaselock loop of LO laser 



A small amount of fundamental light at 1064 nm from the Diabolo laser was split off and interfered with light from 
a Lightwave 126 NPRO 1064 nm laser. The beat note between these two lasers was detected and electronically mixed 
with a tunable local oscillator to create an error signal that was subsequently amplified, filtered and sent to a piezo 
on the lasing crystal of the Lightwave laser, which was capable of tuning the frequency of this laser. By adjusting 
the local oscillator frequency, we were able to phase-lock the two IR lasers to any desired frequency difference. To 
measure cluster states, we locked the two IR lasers at degeneracy, while sending the remaining portion of the Lightwave 
beam through a fiber electro-optic modulator (EOM), creating tunable, high-frequency (< 14GHz) sidebands set at 
ft = Uo ^ {n — 1/2) A, where n=l,...,15 is an integer. This modulated field was then sent through another EOM, 
composed of a RbTiAs04 (RTA) crystal, which applied FM sidebands at 9.2 MHz . This doubly modulated beam was 
split into two, with a small portion being sent to the OPO cavity, where the 9.2 MHz sidebands were used to lock the 
OPO on the first-order high frequency sidebands with a variant of the Pound-Drever-Hall method [ ]. The remainder 
of the beam was sent through a filter cavity of the same FSR as the OPO, which, like the OPO, was locked so as to 
transmit solely the first-order sidebands from the fiber EOM. These sidebands were then used as local oscillator (LO) 
fields in the two-tone balanced homodyne detection (BHD) of each given frequency pair of the OPO. The signal beams 
emitted by the OPO were aligned with the local oscillators by using an alignment beam which was mode-matched 
to the OPO through the output coupler. The reflection of this beam off of the OPO was then interfered with the 
output from the filter cavity. The filter cavity's output was aligned to provide maximum homodyne visibility with 
the alignment beam. A visibility of 97% was obtained for the y polarized beam, while 96% was obtained for the z 
beam. The pump beam was then aligned into the OPO in a way that yielded maximum amplification of the resonant 
modes of the alignment beam. After the pump beam alignment, the alignment beam was blocked to make way for 
the OPO locking beam. 



DETECTION PHASES OF QUADRATURE OBSERVABLES 

Nullifier phases 

Our detection system consisted of two balanced homodyne detectors (BHD), one for the y polarized beams and one 
for the z polarized ones, with all photodiodes measuring both sideband frequencies simultaneously. These nullifiers 
were experimentally measured by BHD with the phase-modulated local oscillator beams transmitted through our 
filter cavity. The fiber EOM modulated the input beam e*'^"^ to output FM beams of the form: 

t-iOo] (14) 



am exp 



t ^ iOo> - am exp 
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where Oq is the phase of the microwave modulation signal at sideband frequency shift (n — 1/2) A. By using a piezo- 
electric transducer (PZT), we then added a phase shift Olo on both sidebands. After reflecting off the PZT mirror, 
the local oscillator beams became 



(15) 



From this, we obtain that the detection phases at w±„ were, respectively. 



9+ = 00 + Olo 

9- = -Oo + 9lo + TT. 



(16) 
(17) 



We now consider the variances of the observables measured by our setup 



A^{g,,gy',e_,e+) =gy [A_nA^_) ^ A^^A^^)] ± + , (18) 

where gy^z are the RF attenuation coefficients applied to the BHD signals of respective polarizations. This is to be 
compared to the nullifiers that must be measured 

A±{0) = A-r,,y{0) T An,y{-e) ± [A.nA^) T AnA-0)] ■ (19) 

It follows that, in order to measure the nullifiers (19), we first needed to realize ^+ = —0- + 2p7r = ^, for A-ri{0) + 
An{—0)^ and 6>+ = —6- + {2p + l)7r = ^, for A_^(^) — An{—6), at each polarization. This leads, respectively, to 
O^^Q = {2p — l)7r/2 to observe A_^{0) and 0^j~q = pir to observe A_{0). As Olo is scanned in the experiment, we 
therefore expected to observe each nullifier to be squeezed alternatively, every 7r/2 rad, which we did, see Fig. 2 in 
main text. Also, do note that ^ Olo- Indeed, the aforementioned particular values of Olo must be plugged into 
Eqs. (16,17) and the FM phase Oo must then be changed by 7r/2 in order to yield the noncommuting quadratures — 
Eqs. (2-5) in main text. We changed Oq by changing the length of the RF cable carrying the FM signal to the EOM. 
However, it is essential to note that the nullifier squeezing (19) is independent of 6>, and the variation of Oq confirmed 
this by leaving the squeezing signal unchanged at all times. 

Secondly, we needed to electronically apply the proper attenuation factors on the two BHD results. The AC output 
signal of each detector was therefore split into two, so that we could make synchronous measurements of the following: 
A±{gz,0;0) and A±{0^ gy;0). These were useful for van Loock-Furusawa criterion measurements as well. In order to 
measure nullifiers, we should therefore set gz = 9y = = 0.618, or -2.1 dB, in those synchronous measurements. 
In practice, we had to take into account a 28% LO loss from the uncoated back surface of a polarizing beam splitter, 
which affected the z detector only. (The quantum signal was not affected by this loss as it was linearized polarized and 
incident at the Brewster angle.) We compensated for this LO power difference by adjusting the attenuation factors 
above by 10 log(l/0.72) = 1.4 dB. Thus, instead of attenuating the z channel by -2.1 dB, it was only attenuated by 
-0.7 dB. Similarly, the y detector was attenuated by -3.5 dB. 



Additional control of nullifiers via polarization 



We introduced an experimental control of the relative phase between the y and z polarized LO beams, by use of 
a half waveplate (HWP) and a quarter waveplate (QWP) in succession. It is possible to ensure that a subsequent 
polarizing beam splitter (PBS) will equally split our LO beam by setting the axis of the QWP to 45°, with respect 
to the vertical. The HWP can then be adjusted to give the desired relative phase between the two LO beams. The 
application of these waveplates on vertically polarized LO beams, with frequencies uJri-,^-n can be explained with 
Jones calculus. We can write the effect of the QWP with axis of transmission at 45° with the vertical axis as Wqwp 
and the effect of HWP with axis of transmission at Ohw p with the vertical axis as Whw p • 

Tj/ fn\ fcos20 sin26> \ 
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Therefore, after effecting Jones matrices on vertical polarized input beam we can write the output beam as: 



1 -A /cos 2(9 sin 2(9 \ /I 
-i 1 j Uin2(9 -cos2(9j lo 



cos 20 — i sin 20 
-i cos 20 + sin 20. 



(22) 



where = Ohwp- Eq. (22) shows that for Ohwp = tt/S, Ohwp = 37r/8, the output beams can be respectively 

written in matrix form as (1 — i) and —(1 -\- i) these particular values of Ohwp^ the local oscillator 

beam is linearly polarized. Changing Ohwp by 7r/4 can apply a tt phase shift between the vertical and horizontal 
polarizations. 



Making simultaneous measurements of the two nullifiers allows us to quickly ensure that Ohw p is set correctly. As 
mentioned before, our homodyne measurement yields nullifiers only for Ohwp = 7r/8+p7r/2 If the angle of the HWP is 
set incorrectly, say with Ohwp = 7r/4, then there is a 7r/2 phase shift between the vertical and horizontal polarizations, 
and the two BHDs will measure quantum amplitudes in quadrature from one another, e.g. A-n.yifi) — ^+n,y{—^) and 
A_^^2(6>+7r/2) — 74+^^2(— ^— 7r/2). These expressions can never correspond to nullifiers, but since we are experimentally 
obtaining finite squeezing, they can, in some cases, give us squeezing which is comparable to the squeezing of the 
nullifiers, and it becomes necessary to distinguish such cases from the nullifiers. In Fig. 7, we plot the measured 
variances in this case. Comparing Fig. 3 in the main text and Fig. 7, it is readily apparent that there is a large 
difference in the relative phases of the two nullifers for ^1,0 = tt/S vs. 6>lo = 7r/4. It is important to remark here 




2.0 3.0 4.0 5.0 
Local Oscillator Phase (rad) 



FIG. 7. Theoretical and experimental variances for Ohwp = 7r/4. The theoretical trace is for r=0.275. 



that, as the amount of squeezing increases, the width of the squeezing hole narrows as the destructive interference 
that leads to the cancellation of the antisqueezed terms in the nullifier variance becomes more stringent. Therefore, 
the aforementioned non-optimum cases, already clearly distinguishable in our present conditions, will become all the 
more obvious. 
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CROSSCORRELATIONS BETWEEN DIFFERENT CLUSTER STATES 



In this subsection, we give the calculation result of the variance of nullifiers (19) that involve modes in two different 
clusters, for example at frequencies ui and uj2. The results are 

AA.{0) = {$ [A^,y{0) + A2,y{-9)] - [A,A9) + A^^-O)]} = + ^) ^2 cosh (^|) (23) 

AA+{0) = A {[A^je) - A2,y{-0)] + * [Ai,,{0) - A2A-m = (1 + *)\/2cosh(2r$) (24) 

As is clear from Eqs. (23,24), writing nullifiers with the wrong modes can never yield squeezing. Figure 4 in the main 
text displays the experimental results which confirm these predictions. We obtained these results by phase- locking 
the two IR lasers at appropriate frequency differences. 



VAN LOOCK-FURUSAWA INSEPARABILITY CRITERIA 



In the situation where there is statistical mixing, we can nonetheless prove complete quadripartite nonseparability 
by showing that our squeezing levels violate the van Loock-Furusawa (vLF) inequalities [6], constructed here from 
variances of A±{gy^gz; —0^ 0)^ where gz and gy denote subtraction gains. In the limit of infinite squeezing, the optimal 
subtraction gains are = = as these coincide with the nullifiers. It is important, though, to point out that 
inseparability is a weaker requirement than state preparation for pure states and the choice of observables above is 
just an optimal one to ensure inseparability for a mixture; it doesn't constitute a rigorous proof that we have created 
a cluster state, only a quadripartite entangled state that should nonetheless be very close to a cluster state given that 
we use cluster-state nullifiers to violate the vLF inequalities. As mentioned previously, the value of phase 9 Eq. (19) 
does not change the squeezing. 

In order to derive the sufficient conditions for full inseparability, we may just consider all possible bipartitions of 
the four modes, and for each case, we enumerate the necessary conditions for separability. As the contraposition, we 
obtain the sufficient condition for full inseparability. 



One mode - three mode separability 



According to the vLF criteria, if the first mode is separable from the other three modes, the following inequalities 
are satisfied, 

{Q-{g,f) + {P-{g.f) > ^[|l X 1| + |1 X (-1) + (-5.) X (-5.) + {-g,) x g,\] 

= 1, (25) 
{Q+{gy)^) + {P+{gyf) ^ g% (26) 
{Q-{g.f) + {P+{gy?) (= {Q+{gy?) + {P-{g.?)) > gy, (27) 

where Q-{gz) = A_(0, 0, 0), etc. The former two of these three equations can be simplified as 

2{Q_{g,f) ^ 1, (Al) 

2{Q+igyf) > gl (A2) 

because of the phase insensitivity of nullifiers, namely (uf) = (vf) with i = {1,2}. Equation (27) can be generalized 
as {Q-{gz)'^) X a + (P+(^^)^)/a ^ gy with any real number a, or equivalently. 



2^{Q.{g,Y){Q+{gyY) ^ gy. (A3) 
Similarly, if the third mode is separable from the other three modes, the following are satisfied, 

2{Q.{g,f)^gl (Bl) 

2{Q+{gyf) ^ 1, (B2) 



2J{Q-{g.)^){Q+{gyf)^g.. (B3) 
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In the case that the second (fourth) mode is separable, the same inequahties for separabihty of the first (third) mode 
are satisfied. 



Two mode - two mode separability 

Next, we divide the four modes into two sets of two. If the first and second modes are separable from the other two 
modes, 



2y/{Q.{g,y){Q+{gyr)^g,+gy. (C) 

Here, we neglect the other two inequalities because they are trivial (> 0). In the case that the first and third modes 
are separable from the other two. 



2{Q_{g,f)^l + gl, (Dl) 



2{Q+{gyf)^l+gl (D2) 



2^(Q_(fl,)2)(Q+(5,)2) ^\g,-gy\. (D3) 
In the case that the first and fourth modes are separable from the other two, 

2{Q_{g,f)^l-gl (El) 

2{Q+{gyf)^l-gl (E2) 



2^(Q_(g,)2)(Q+(5,)2) ^ - (E3) 

Sufficient conditions for inseparability 

Sufficient conditions for full inseparability can then be obtained by violating one inequality from each group. The 

necessary squeezing levels for inseparabilty can be found by normalizing each of the above inequalities by the shot 
noise level of This yields: 

{Q-{9.f)N < (Al') 

{Q+{gyf)M < (A2') 



J{Q.{9,f){Q+{gyf) < (A3') 
^ ^ ^(l + gl){l+gl) 



9l 



{Q-{9.r)N < jf-^, (Bl') 
{Q+{gy?)N < (B2') 



J{Q.{g,f){Q+{9yY) < (B3') 

" ^ ^d + gD^l+gl) 

J{Q-{g.Y){Q+{gyY)^ < . (C) 
^ ^ ^d + gDil+gl) 
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{Q-{g,f)N < 1, (Dl') 



{Q+{gyf)N < 1, (D2') 



{Q_(5,)2)(Q+(5,)2) < , (D3') 

^ y'(l+52)(i+52) 



(Q_(5,)2)^ < i-g, (El') 
{Q+{gyf)N < (E2') 



y 



{Q.{9.f){Q+{gyY)^^ < , _ , (E3') 



For ^2 = = 1/^5 where most of our data was taken, it is necessary that we have 



{Q.{g,r)N{Q+{gyr)N < ^ (28) 



and any of the fohowing 



^2^ , a/5-1 



{Q-{g.r)N < (29) 
{Q^{gyf)N < (30) 



'{Q-{g.)')N 

{Q+{9y?)N < 



(31) 
(32) 



where the subscript N denotes normahzation with the shot noise level {^f^ = ^^). Equation (28) requires that 
our two nullifers have a mean squeezing value of -0.5 dB, which we easily exceed. Meanwhile, equation (32) requires 
-1.4 dB and -3.5 dB squeezing relative to the shot noise level, which we have experimentally achieved for all measured 
sidebands, thereby ensuring quadripartite entanglement. We have enough squeezing to easily satisfy the inequalities 
given in (B2'), (C), (Dl'), and (D2'). The violations of (Al') and (E2') are, however, much smaller. It is possible to 
observe larger violations of the van Loock-Furusawa criteria in these cases by adjusting the gains away from 1/$. By 
setting gz = 9y = 0-3, we reduce the necessary squeezing for (Al') from -1.4 dB to -0.4 dB. Similarly, the necessary 
squeezing for (E2') is reduced from -3.5 dB to -0.8 dB. As can be seen in Fig. 6 of the main paper, the squeezing 
levels at these gains are largely exceeded. 
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